Non-linear coupling between the two oscillation modes of a dc-SQUID 
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We make a detailed theoretical description of the two-dimensional nature of a dc-SQUID, analyz- 
ing the coupling between its two orthogonal phase oscillation modes. While it has been shown that 
the mode defined as "longitudinal" can be initialized, manipulated and measured, so as to encode a 
quantum bit of information, the mode defined as "transverse" is usually repelled at high frequency 
and does not interfere in the dynamics. We show that, using typical parameters of existing devices, 
the transverse mode energy can be made of the order of the longitudinal one. In this regime, we can 
observe a strong coupling between these modes, described by an Hamiltonian providing a wide range 
of interesting effects, such as conditional quantum operations and entanglement. This coupling also 
creates an atomic-like structure for the combined two mode states, with a V-like scheme. 



PACS numbers: 85.25.Cp, 03.67.Lx 

The control of physical systems at the quantum level 
represents a significant technological challenge. But, this 
goal is of extreme interest since successful applications 
would allow us to probe fundamental aspects of quan- 
tum physics and to demonstrate the claimed advantages 
of quantum computation. A number of physical sys- 
tems have shown, in the last decade, their potentiality 
as "quantum devices" . Examples are atoms in a high fi- 
nesse cavity [H , trapped ions Q and solid state quantum 
circuits [^^. 

Superconducting circuits (SC) appear as suitable can- 
didates to study controllable macroscopic quantum sys- 
tems. During the last decade, many experiments on such 
devices have demonstrated effects previously observed on 
atomic and quantum optics experiments [3 0] . Further- 
more, the extreme flexibility in choosing the characteris- 
tic parameters of SC have made possible the realization of 
new interaction regimes, such as the ultra strong coupling 
regime 0, Si- these experiments, SC act as artificial 
atoms with one degree of freedom, described either as a 
2-levcl system or a multilevel system. 

In this work, we show that a dc-SQUID with a large 
loop inductance is an artificial atom with two degrees of 
freedom. We make a theoretical study of the two or- 
thogonal phase oscillation modes and analyse the cou- 
pling between them. One of the modes, called hereafter 
the longitudinal mode (LM), was extensively studied and 
demonstrated a large variety of quantum phenomena . 
For instance, using the two lowest energy states, all re- 
quired operations necessary to define a qubit (coherent 
manipulation, initialization, coupling and measurement) 
have been performed [lo| . The quantum dynamics of the 
transverse mode (TM), which is perpendicular to LM, 
have never been studied theoretically or experimentally, 
although it presents appealing potential advantages, such 



as longer coherence times and an optimal point for energy 
relaxation. 

Moreover, we also show that the nature and the 
strength of the coupling terms between the two oscillation 
modes can be engineered to a large extent. They lead 
to a wealth of novel effects, that were initially demon- 
strated or proposed in atomic physics. Examples include 
quantum non-demolition measurements, the realization 
of conditional quantum operations, the controlled cre- 
ation and measurement of entanglement between the two 
mode. In addition, the coupling term in the dc-SQUID 
can result in a V-type 3-level structure, that has recently 
attracted a lot of interest . All these potential exper- 
imental achievements involve realistic circuit parameters 
and should be observable with present day technology. 

A dc-SQUID is composed of two Josephson junctions 
(J J) with critical current /q, capacitance Cq and phase 
difference (t>i and 02 , embedded in a loop of inductance 
L (Fig. [T^). The phase dynamics of this system can 
be mapped on the evolution of a fictitious particle of 
mass m = 2Co (^o/Stt), with two degrees of freedom 
X = (01 -I- 02)/2 a nd y — (0i — 4>2)/'^, which is expe- 
riencing a potential (Hill: 

U(x,y) = C/o[- cosxcosy - —x + b[y - -— ], (1) 

Zlo V 5>o / 

where (f>o = h/2e is the superconducting fiux quantum, 
Uq = 2/o ($o/27r), lb and <!>{, are respectively the cur- 
rent and flux bias and b = the ratio between 
the Josephson inductance and the loop inductance. For 
Ih < 2/o, this potential presents a periodic series of 
equivalent local minima separated by potential barriers 
(Fig.[T]D). For simplicity, we consider a perfectly symmet- 
ric dc-SQUID, with symmetric inductances, JJs critical 
currents and capacitances. 
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X// (longitudinal mode) 



Yx(transverse mode) 



FIG. 1: (a) Electrical scheme of an inductive dc-SQUID 
(b) 2D potential for a dc-SQUID with b = 3, h ^ and 
$6 — 0.51$o; M and M' indicate two diflerent local minima 
separated by a saddle point C. The dashed black line repre- 
sents the trajectory of minimum energy, (c) potential of the 
LM (left side) and TM (right side), related to the motion of 
the fictitious particle trapped at the local minimum M. 



To describe the 2D dynamics of the particle trapped 
in a potential minimum, we introduce the basis (x||,xj^) 
corresponding to the motion parallel (||) or transverse 
(_L) to the trajectory of minimum energy [l^ . It is ob- 
tained from a rotation with an angle 9{Ib, ^b) of the ba- 
sis (x, y). The system dynamics is then governed by 
the Hamiltonian "H = "Hy + C + TL± which describes 
two anharmonic oscillators whose interaction is medi- 
ated by a coupling term C (Fig. [TJ;). Up to order 4, 
Tia = huJa X [^{pI^ + xl) - aaX^ + Sax'^] with a =11 or ±. 
Here, Xa and pa are the reduced (dimensionless) position 
and momentum operators in the direction a; they satisfy 
the commutation relation [xq.,Pq] — il; the oJq. are the 
characteristic frequencies of each mode; in general, they 
are different with wy < ujj_. For relevant bias points, the 
oscillators are weakly anharmonic: Sa ^ 1 and CTq, ^ 1. 
Thus the energy En^ associated with the eigenstate \na) 
of Ha can be obtained from a simple perturbative ap- 
proach. Specifically, — £'ri„-i = {l — Aana)hLUa 
with Aa = ^cTa +2>5a [10| . The coupling Hamiltonian 



reads C = fi(?2iiya;_L + fi'gi2X\\x\ -\- hg22x'^\x\. Note that 
here the linear coupling term xyij^ is cancelled by the 
(x||,x_l) basis choice. The fourth order coupling terms 
x\\x\ and xj|i_L are negligible at finite bias current and go 
to zero at zero bias current. We have also the frequency 
hierarchy 321,^12,522 < 

The LM alone was previously studied [13, lil- 17|. 
Quantum state detection was realized using a nanosecond 
pulse to produce a state dependent escape rate, through 
the potential barrier. At zero-current bias, the LM mode 
is insensitive to current noise Nevertheless, away 

from the zero- flux bias point, it is still flux noise sensitive 
as other phase qubits [l8|. Until now, the TM has not 
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FIG. 2: Energy spectrum. Transition frequencies between 
ground state |0||,0x) and excited states |n||,nx) versus 
at zero current bias of the artificial atom described by two 
oscillators without coupling (dashed lines) and taking into 
account the coupling (continuous lines) with /o = 710 uA, 
Co = 510 fF, 6 = 0.73. The two solid blue curves highlight 
the states |0||, 1±) and |2||,0±), chosen to illustrate the effect 

of the coupling term n^f^ . The area of maximum entangle- 
ment between these two states is found inside the dotted gray 
rectangle. 



been studied and the resonant dynamics of the dc-SQUID 
have been restricted to one dimension only. Indeed the 
plasma frequency of the TM strongly increases with h. 
In previous dc-SQUID studies, the loop inductance was 
significantly smaller than the Josephson inductance, i.e 
6 > 1, leading to a high TM plasma frequency [13, [isj . 
By making the loop inductance larger than the Joseph- 
son inductance one can decrease the TM frequency and 
make it closer to the LM one, as shown in Fig. [2l The 
TM presents several interesting theoretical features that 
demand to be experimentally investigated and exploited. 
For instance, current fluctuations induced by the envi- 
ronment lead to a TM relaxation rate F^ = gi{d)Si{uj±) 
where Sj {uj± ) is the spectral density of the current noise 
and gi{0) the coupling strength between current and the 
TM. For a bias current /(, = 0, this coupling goes to 
zero for the TM. Thus, a longer relaxation time is ex- 
pected. Moreover, the dephasing time depends on ^f^, 
which goes to zero at = 0, and on which is one 
order of magnitude smaller than for the LM (see Fig.[2|). 
Therefore, the TM is expected to suffer from a smaller 
dephasing than the LM. 

In addition, the coupling between the TM and the 
LM leads to a number of physical effects, expanding the 
range of applicability of SC devices in quantum informa- 
tion field. Hereafter we will concentrate our study on 
these coupling effects. The Hamiltonian of the uncou- 
pled degrees of freedom is given by Ho — Hn + 'H±. 
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The energies E„ 



E„ . + En^ of its eigenstates 



= Ww) ® |n_L) are shown in Fig. 2 as dashed 
hnes, for a particular set of circuit parameters. Since 
522) ^12) <72i ^ '^±7^11, the interaction term C is treated 
as a perturbation. To perform a detailed analysis, we 
introduce the annihilation operators a and b of the LM 
and TM, and the occupation number operators ny = a^a 
and nj^ = h^h. One then obtains C = 2^2 ^■'^^^^^^ ^ 
+ 2fi|| + 1) X (6 + 6t) + ^hgi2{y^ + + 2n^ + 

1) X (a + at) - \hg22{a^ + + 27i|| + 1) x + 6^^ + 
2n±^ + 1). In the interaction picture, the interaction 
Hamiltonian is Hint = UCU\ with U e-^('^ii+«^)*. 
The effect of the coupling is captured by the matrix 
elements , n^|7^i,it l^^n , 'T'i) = 7 
where he = E,-,-^^ - Er^ + -En'^ - En^. 

In order to couple both modes in a non negligible 
way, {n'^<^Tn'^_\'Hint\n\\,n±) must oscillate at a frequency 
e much slower than the coupling strength between the 
two states. Only two terms can potentially verify this 



condition in C: % 



and % 



(1) 



eff - 3%22 (2ni| + 1) x {2n^ + 1) 
J§ ^ i^92i{ba^'^ + b^a^)- The first term is 
present independently of the values of riy and nj^. It 
induces a conditional energy shift for both modes. The 
second term leads to a single desexcitation of the TM 
and the simultaneous double excitation of the LM (and 
the hermitian conjugated process). It is relevant only 
when a quasi-resonant condition is verified between states 
|7T,||,7ij^) and |n|| =p2,nj^ ± 1). The deviation from the 
resonant condition between these two states is quanti- 
fied by e„||,„^($6) = (£^„,| + E^^ - En^^-2 - En^+i)/h w 
2a;j| — — (2n|| + 1)A||W|| + nj_Aj_CL'j_. One can choose 
a set of circuit parameters so that 2lo\\ k, u)_\_. As shown 
in Fig. [21 tuning W|| and a;_L with the bias flux, we cither 
satisfy the exact resonant condition {en^^,n±{^b) = 0), 
or ensure an off-resonant condition (leny .n^ ('I?;))! ^ 52i)- 
We now study the rich physics associated with "Hgg and 

(2) 

"Hgff . Depending on the experimental parameters, each 
of these two terms can be chosen to dominate the cou- 
pling Hamiltonian. 

We start with Hg^, setting \en^^,n_i_{^b)\ ^ 321, so that 

V}-^^ can be neglected. H^^^ is very similar to the Kerr ef- 
fect Hamiltonian [l9| . with typical values of g22/27r rang- 
ing from 50MHz to 500MHz. It allows i) the realization of 
quantum non-demolition (QND) measurement of the LM 
or TM [iOl, ii) conditional quantum logic (CQL) gates 
2l| and iii) the realization of a V-type 3-level atom, in- 
troduced initially in atomic-CQED. For example, QND 
measurement of the LM can be realized by TM spec- 
troscopy: 'm}'^^ shifts in frequency the transition between 
the Ho eigenstates |n||,0_L) and |n||,l_L) proportionally 
to (2?i|| + 1). Using the two lowest energy states of both 
modes as qubits, this frequency shift can also be used 
to implement CQL gates. A microwave 7r-pulse resonant 
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FIG. 3: (color online) Population transfer between states of 
the TM and the LM as a function of time for A = 0.69(?2i- 
Time is in units of gj"^^ . The blue continuous line shows the 
population of state |0||, Ix), set equal to 1 at t = 0, and the 
pink dashed line is the population of state |2||,0x), initially 
equal to 0. The inset shows the absolute value of the negativ- 
ity as a function of time, providing evidence that ma:ximally 
entangled states between the two modes can be created. 



with the |0||,lx) ^ longitudinal transition can 

switch both states without affecting |0||,0x) and |l||,Ox)- 
Therefore the qubit defined by the TM acts as a control 
qubit in a CNOT gate. Such gate have been successfully 
realized with fiux qubit circuit (22j , but away from the op- 
timal point. In the proposed dc-SQUID, such operations 
can be realized at the optimal point (/^ = O,^}, = 0), in- 
sensitive to flux and current noise in first order. Finally, 
the three states |0||,0j^), |l||,Oj^) and |0||,lx) reahze a 
V-like level structure, that was also proposed recently 
in charge qubit [2F 
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inserting 



As proposed in Ref. 
such a V-type atom in a microwave resonator opens av- 
enues for a wealth of effects discussed up to now in the 
context of atomic physics, such as quantum beats [23|, 
the realization of dark states and the irnplementatiou of 
non-linear quantum field Hamiltonians |28| . 

We now discuss H^.^^ in the case ^22 ^ 521 ■ so that 
"h!"^^ can be neglected. Experimentally, this condition is 
satisfied in a dc-SQUID with the parameters presented 
in Fig. [5] In this case, couples states |n||,nx) to 

states |n|| T2,nj_ ± 1) when |e„|| ($h)| <C .921- In the 

f2) 

following, we show that H),^ can be used to i) create the 
Fock state |2||) in the longitudinal mode, starting from a 
single excitation quanta in the transverse mode, and ii) 
entangle the states of the two oscillation modes. 

To illustrate the first point, we consider the subspace 
spanned by the pair of states |2||,0x) and |0||, The 
dependences of the associated energies on the bias flus are 
shown in Fig.[2]as solid blue curves. These two states sat- 
isfy the resonance condition at a bias flux <I>r = 0.54$0j 
and 'h!'^^ then results in an avoided crossing (inside the 
gray rectangle). Away from the resonance, the coupling 
is negligible and the system eigenstates are those of the 
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uncoupled Hamiltonian Hq. Initially, $b is set to and 
a TT-pulse drives the TM from |0||,0_l) to |0||,1_l) (this 
state has a lower energy than |2||,0j^)). Increasing $b 
slowly enough, it is then possible to realize an adiabatic 
passage through the avoided crossing. This results in the 
conversion of the state |0||,1_l) into the state |2|[,0_l). 
To precise the adiabaticity condition, we suppose that in 
the vicinity of the resonance, the frequency difference be- 
tween the states varies linearly in time: e„||_„^ ($f,) = Xt. 
The adiabatic regime is obtained for a frequency sweep 
parameter A ^ (7r/2)g|j. 

This process is illustrated in Fig. |3] using realistic ex- 
perimental parameters. The initial frequency splitting 
between the two involved states is e2||,0j_(0)/27r « 5 GHz 
and (721 /27r w 600 MHz close to the level anti-crossing. 
With A = 0.69(721, corresponding to a variation of 1.8GHz 
per nanosecond, one can perfectly transfer the population 
from |0||,lj^) to |2|[,0j^). Experimentally, this frequency 
sweeping speed can be realized using a pulse generator 
with a typical rise-time of 5 ns. This time should be 
compared to the coherence damping time Ti. For a sin- 
gle mode dc-SQUID in the state Ti w 100 ns 
Since we expect the TM to present longer values of Ti, 
these relaxation times arc more than enough to perform 
the state transfer protocol. 

Here, we must stress that state transfer is made possi- 
ble by the adiabatic evolution of In the instantaneous 
limit (A ^ (7r/2)g|]^), if the system is initially prepared in 
|0||, 1_l), for instance, it does not have the time to couple 
to state |2|[,0_l) during the fast evolution. Consequently, 
coupling can be disregarded during the flux time vari- 
ation. Such a condition is met using a pulse generator 
with a typical 500 ps rise-time. 

Inspecting again Fig. [3] around the resonance, one ob- 
serves strong oscillations in the population of jOy , 1_l}. As 
we approach resonance, |0|| , Ij^) is no longer an eigenstate 
of Hq. Interestingly enough, the system eigenstates are 
in fact entangled states of the two oscillation modes. The 
inset of Fig. |3] shows the time evolution of entanglement, 
quantified by the negativity TV [l^l- During the adia- 
batic evolution, the system passes through a maximally 
entangled state, evidenced by A/" = 1/2. This considera- 
tion, together with the time evolution regimes explored 
in the previous paragraphs, suggest a procedure to create 
and detect controlled entanglement between the longitu- 
dinal and transverse modes. To prepare a maximally 
entangled state, for instance, one can stop the adiabatic 
evolution at the time corresponding to M = 1/2. By 'in- 
stantaneously' restoring the system to the non-resonant 
regime, it remains a maximally entangled state. One can 
then measure this state and its entanglement by measur- 
ing the population and coherences of the TM and the LM. 
The same type of analysis can be extended to other, more 
excited, TM-LM states, where coupling between modes 
can also occur depending on the value of e„||_„^ ($f,). 

The double excitation of one mode and the simultane- 



ous single desexcitation of another mode (and the hermi- 
tian process), has been experimentally observed in cavity 



QED systems [2j|. In this system, the two modes of a 



high quality microwave cavity (formally equivalent to two 
harmonic oscillators) interact, mediated by an atom, via 
a third order process in perturbation theory. In trapped 
ion systems, such Hamiltonians can also be produced 
by tuning a laser to specific sideband frequencies cor- 
responding to the excitation of the two-dimensional vi- 
brational motion [25j . Usually, the coupling constant is 
proportional to the fourth power of the Lamb-Dicke pa- 
rameter, so the coupling strength is significantly reduced. 
In the present system the non linear coupling is observed 
in the strong coupling limit, i.e. g2i ^ 2ttT^2- 

In conclusion, we have shown that a dc-SQUID is an ar- 
tificial atom with two coupled degrees of freedom, the LM 
and TM. This coupling is fundamentally interesting and 
can be described by a number of effective Hamiltonians 
whose strengths vary according to experimentally tun- 
able parameters. Each Hamiltonian gives rise to a num- 
ber of applications, such as QND measurement, quantum 
logic, entanglement creation and the introduction of a 
more complex level structure to a superconducting arti- 
ficial atom. Further applications arise when we consider 
coupling this system to a microwave cavity, opening the 
way to generalizing a larger collection of results already 
existing in cavity QED. 
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